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Abstract 

The light-front approach is a relativistic quark model and offers many insights to the internal 
structures of the hadronic bound states. In this study, we apply the covariant light-front approach to 
ground-state heavy quarkonium. The pesudoscalar and vector meson decay constants are discussed. 
We present a detailed study of two-photon annihilation P — > 77 and magnetic dipole transition 
V — > -P7 processes. The numerical predictions of the light-front approach are consistent with the 
experimental data and those in other approaches. The relations of the light-front approach with 
the other methods are discussed in brief. 
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I. INTRODUCTION 



Heavy quarkonium provides a unique laboratory to study Quantum Chromodynamics (QCD) 
for the bound states of heavy quark-antiquark system. The nature that the heavy quarkonium is 
relevant to a non-relativistic treatment had been known for a long time [1.]. Although the non- 
relativistic QCD (NRQCD), an effective field theory, is a powerful theoretical tool to separate the 
high energy modes from the low energy contributions, the calculations of the low energy hadronic 
matrix elements rely on model-dependent non-perturbative methods in most cases. Prom the point 
of view of the non-perturbative QCD, there is no one method which is uniquely superior over the 
others. Many methods were employed in heavy quarkonium physics, such as lattice QCD, quark- 
potential model, etc. (for a recent review see [2J). The light-front quark model, in which a hadronic 
matrix element is represented as the overlap of wave functions, offers many insights into the internal 
structures of the bound states. In this study, we will explore the heavy quarkonium from a quark 
model on the light front. 

The light-front QCD has been developed as a promising analytic method for solving the non- 
perturbative problems of hadron physics |3]. The aim of the light-front QCD is to describe the 
hadronic bound states in terms of their fundamental quark and gluon degrees of freedom. It 
may be the only possible method that the low energy quark model and the high energy parton 
model can be reconciled. For the hard processes with large momentum transferred, the light- 
front QCD reduces to perturbative QCD (pQCD) which factorize the physical quantity into a 
convolution of the hard scattering kernel and the distribution amplitudes (or functions). In general, 
the basic ingredient in light-front QCD is the relativistic hadron wave functions which generalize 
the distribution amplitudes (or functions) by including the transverse momentum distributions. It 
contains all information of a hadron from its constituents. The hadronic quantities are represented 
by the overlaps of wave functions and can be derived in principle. 

The light-front quark model is the only relativistic quark model in which a consistent and fully 
relativistic treatment of quark spins and the center-of-mass motion can be carried out Q]. This 
model has many advantages. For example, the light-front wave function is manifestly Lorentz 
invariant as it is expressed in terms of the internal momentum fraction variables which is inde- 
pendent of the total hadron momentum. Moreover, hadron spin can also be correctly constructed 
using the so-called Melosh rotation. This model had been successfully applied to calculate many 
phenomenologically important meson decay constants and hadronic form factors 

On the light front, the non-relativistic nature of a heavy quarkonium is represented by that the 
light-front momentum fractions of the quark and antiquark is close to 1/2 and the relative trans- 
verse and the z— direction momenta are much smaller than the heavy quark mass. The Lorentz 
invariant light-front wave function and the light-front formulations provide a systematic way to 
include the relativistic corrections. There is no conceptual problem to extend the light-front ap- 
proach into the heavy quarkonium. We will apply the covariant light-front approach 0, @] to 
the ground-state s-wave mesons which include 1 5o pseudoscalar mesons (P) rjc, % and 3 Si vector 
mesons (V) J/ift,T(lS) as our first-step study along this direction. The main purposes of this 
study are threefold: (1) Is the light-front approach applicable into the heavy quarkonium? In 
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concept, the light-front quark model is the relativistic generalization of the non-relativistic quark 
model. The phenomenological success of the previous non-relativistic quark-potential model should 
be reproduced in the light-front approach. In particular, we will examine the validity of the light- 
front approach in three types of quantities: decay constants, two-photon annihilation P — > 77 and 
ma gnet ic dipole transition V — > P7. In most literatures, these processes were explored separately 
flCl Jill, [12I 13]. To study them simultaneously can better constrain the phenomenological param- 
eters and check the consistency of the theory predictions. (2) The 775 meson has still not been 
observed in experiment [3]. We will present the numerical prediction for the branching ratios for 
Vb 77 an d T — * 7^7 processes. (3) What is the relation of the light- front approach with the other 
approaches? In the non-relativistic approximations, the light-front approach will be closely related 
with the non-relativistic quark-potential approach. For the process of P — > 77 which is light-front 
dominated, the light-front approach reduce to the model-independent pQCD. 

The paper is organized as follows. In Sec. II, we give a detailed presentation of the covariant 
light-front approach for heavy quarkonium. It contains a brief review of the light-front framework 
and the light-front analysis for the decay constants of P and V mesons and the processes P — > 77, 
V — * Pj. In Sec. Ill, the relations of the light-front approach with the non-relativistic approach 
and pQCD are discussed. In Sec. IV, the numerical results and discussions are presented. Finally, 
the conclusions are given in Sec. V. 



II. FORMALISM OF COVARIANT LIGHT-FRONT APPROACH 
A. General formalism 

A heavy quarkonium is the hadronic bound state of heavy quark and antiquark. In this system, 
the valence quarks have equal masses mi = m 2 = rn with m the mass of heavy quark c or b. Thus 
the formulae associated with the term (mi — m 2 ) vanish and will lead to some simplifications. In 
this section, we will give the formulae specially for the quarkonium system. 

The momentum of a particle is given in terms of light-front component by k = (k~,k + ,k±) 
where k^ = k° db A; 3 and k± = (k 1 , k 2 ), and the light-front vector is written as k = (k + , k±). The 
longitudinal component k + is restricted to be positive, i.e., k + > for the massive particle. By 
this way, the physical vacuum of light-front QCD is trivial except the zero longitudinal momentum 
modes (zero- mode). We will study a meson with total momentum P and two constituents, quark 
and antiquark whose momenta are p\ and p 2 , respectively. In order to describe the internal motion 
of the constituents, it is crucial to introduce the intrinsic variables (xi,p±) through 

pf = x 1 P + , p 1± = x 1 P± + p±; 

p+ = x 2 P + , P2± = x 2 P±-p±, (2.1) 

where Xi are the light-front momentum fractions and they satisfy < xi,x 2 < 1 and x\ + x 2 = 1. 
The invariant mass Mq of the constituents and the relative momentum in z direction p z can be 
written as 

2 = pi + m^ = _ 1 

xix 2 2 
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The invariant mass Mq of qq is in general different from the mass M of meson which satisfies M 2 = 
P 2 . This is due to the fact that the meson, quark and antiquark can not be on-shell simultaneously. 
The momenta p± and p z constitute a momentum vector p = (p±,p z ) which represents the relative 
momenta in the transverse and z directions, respectively. The energy of the quark and antiquark 
ei = 62 = e can be obtained from their relative momenta, 



lm 2 +p\+p 2 . (2.3) 
It is straightforward to find that 

e-p z e + p z M , . 

x\ = , xo = , e = . 2.4) 

2e 2e 2 y J 

To calculate the decay constants or decay amplitude, the Feynman rules for the vertices of 

quark-antiquark coupling to the meson bound state are required. In the following formulations, we 

will follow the notations in Q]. The vertices Tm for the incoming meson M are given as 



Hpjb for P; 

7,4 - T7j-{pi ~ P2)y, for V. (2.5) 



Wy 

After performing a one-loop contour integral to be discussed below which amounts to make one 
quark or antiquark on its mass-shell, the function Hm and the parameter Wy are reduced to Km 
and wy, respectively, and they are written by 

h P = hy = (M 2 - Mo) /^ J M ^X2,P±), (2.6) 

and 

wy=M + 2m. (2.7) 

The form of the function Km and the Feynman rule for Tm are derived from the light-front wave 
function which describes a meson bound state in terms of a quark q\ and an antiquark qi- The 
light-front wave function contains two parts: one is the momentum distribution amplitude (fi(x2,p±) 
which is the central ingredient in light-front QCD, the other is a spin wave function which constructs 
a state of definite spin (S, S z ) out of light front helicity eigenstates (Ai, A2). The spin wave function 
is constructed by using the Melosh transformation and its spin structure has been contained in Eq. 
(|231) . 

The momentum distribution amplitude 4>(x2,p±) is the generalization of the distribution am- 
plitude (f)(x) of the pQCD method and can be chosen to be normalizable, i.e., it satisfies 

dxd 2 p± 2 



2(2?r) 3 A^)\' = 1- (2-8) 

In principle, 4>{x2,Pa_) is obtained by solving the light-front QCD bound state equations -ff£ir|\E f ) = 
M\^>) which is the familiar Schrodinger equation in ordinary quantum mechanics and Hlf is the 
light-front Hamiltonian. To see the explicit form of the light-front bound state equation, let us 
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consider a quarkonium wave function. The light-front bound state equation can be expressed as: 



M 1 



E 



( \ 



my 



/ (qq\H int \qq) 
(qqg\H int \qq) 

V '■ 



V i / 

(qq\H int \qqg) 
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(2.9) 




Of course, to exactly solve the above equation for the whole Fock space is still im pos sible. Currently, 
two approaches have been developed. One is given by Brodsky and Pauli 
called discretize light-front approach, the other by Perry, Harindranath and Wilson 
based on the old idea of the Tamm-Dancoff approach 22J, [23J that truncates the Fock space to only 
include these Fock states with a small number of particles. Furthermore, if one can eliminate all 
the high order Fock space sectors (approximately) by an effective two-body interaction kernel, the 
light-front bound state equation is reduced to the light-front Bethe-Salpeter equation: 

dyd 2 k\ 



M 



k\ + m 2 
x(l — x) 



-V ef f(x, k±,y, k' ± )y g g(y, k' ± 



(2.10) 



2(2tt) 3 

One may solve the Bethe-Salpeter equation for finding the relativistic bound states. However, the 
Bethe-Salpeter equation only provides the amplitude of a Fock sector in the bound states so that 
it cannot be normalized. In other words, the Bethe-Salpeter amplitudes do not have the precise 
meaning of wave functions for particles. In addition, the advantage for Eq. (|2.9p with the Tamm- 
Dancoff approximation is that it provides a reliable way to study the contribution of Fock states 
which contain more particles step by step by increasing the size of truncated Fock space, while the 
Bethe-Salpeter equation Eq. (|2.10p lacks such an ability. Some studies on nonperturbative features 
of light-front dynamics were focused on the 1+1 field theory. Typical examples are: the discretized 
light-front quantization approach for the bound states in the 1+1 field theory developed by Pauli 
and Brodsky 15|, [2J] , the light-front Tamm-Dancoff approach for bound state Fock space truncation 



discussed by Perry et al. [19j, |20j]. However, at the present time, how to solve for the bound states 
from 3+1 QCD is still unknown. We are satisfied with utilizing some phenomenological momentum 
distribution amplitudes which have been constructed phenomenologically in describing hadrons. 
One widely used form is the Gaussian-type which we will employ in the application of covariant 
light-front approach. 



B. Decay constants 

In general, the decay constants of mesons fpy are defined by the matrix elements for P and V 
mesons 

<0|^|P(P)> = i/ P P M , 

(0\V„\V(P)) = M v f v e^. (2.11) 
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FIG. 1: Feynman diagram for meson decay constants, where P is the momentum of meson, p\ is the 
quark momentum, p 2 is the antiquark momentum and T denotes the corresponding V-A current. 



where P^ is the momentum of meson and e M is the polarization vector of V meson. The Feynman 
diagram which contributes to fpy is depicted in Fig. [TJ The meson decay constant plays an 
important role in determining the parameters of the distribution function 4>(x2,p±), in particular, 
the quark mass and a parameter (3 characterizing the hadronic "size" for a Gaussian wave function. 
The decay constants have been calculated in 0, [S] and are the same as our results. Thus we simply 
provide the formulae for fpy here. 

For a pseudoscalar quarkonium, the decay constant is represented by 



fF = ^ J dX2d ^^M MX ^ P±) 



v^Vc 



8vr 3 



m 



dx 2 d p± =4> P (x 2 ,p±). (2.12) 
/ m 2 + p\ 

where iV c = 3 is the color number and m denotes the mass of the heavy quark. In Eq. (|2.12p . we 
have used the relation 

Mov^I^ = \/m 2 +pl (2.13) 

for a quarkonium . 

For the vector meson, the decay constant in the covariant approach is represented by 

4>v(x2,P±)- (2.14) 



V2N C f 2 
fv= 8^Mj dX2d P± 



m 2 +p 2 j_ 



Xl M 2 -pl + —p\ 



Eq. (|2.14p coincides with the result in [7j when m\ = mi- Note that the 1/wy part of Eq. 
(|2.14p is different from that in the conventional approach, for example, 0]. The reason is that the 
conventional approach is not covariant and contains a spurious dependence on the orientation of the 
light front. The relevant calculations are not free of spurious contributions for transitions involving 
vector meson. Zero modes, which relate to the p~ integration for p + = 0, are required to eliminate 
the spurious dependence and contribute to the 1/wy part of Eq. (|2.14p . More detailed discussions 
about this point can be found in 0, B|. The decay constant fy is related to the electromagnetic 

47r a 2 2 




decay of vector meson V — > e + e 
where cy is factor related to the electric charge of the quark that make up the vector meson. 



) = TS h»#. (2.15) 
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C. P -» 77 



Charge conservation requires charge conjugation C = +1 state coupling to two photons. Thus 
only pseudoscalar meson can transform into two photons while the vector meson is forbidden. In 
the process of P — > 77, the final two photons are both on-shell. For the purpose of illustration, it 
is useful to consider a more general process P — ► 77* with one photon off-shell. We introduce a 
transition form factor Fp 1 {q 2 ) arising from the -P77* vertex. The P — > 77 process is related to the 
form factor at q 2 = 0, i.e., i ? p 7 (0). The form factor Fp 1 (q 2 ) is defined by 

Afi = -i e 2 F P ^q 2 )e ilvpa P v q{e (J . (2.16) 

where A^ is the decay amplitude of the process P — * 77* and 91(e) the momentum (polarization) 
of the on-shell photon. 



7*(g) 




(a) (b) 

FIG. 2: Feynman diagram for P — > 77* process where P in the parenthesis denotes the momentum 
of meson. The diagram (b) is related to (a) by the exchange of two photons. 



The transition amplitude for the process of P —* 77* can be derived from the common Feynman 
rules and the vertices for the meson-quark-antiquark coupling given in Eq. (12. 5|) . In the covariant 
light-front approach, the meson is on-shell while the constituent quarks are off-shell and the mo- 
mentum satisfies P = pi + P2- To the lowest order approximation, P — > 77* is a one- loop diagram 
and depicted in Fig. [2j The amplitude is given as a momentum integral 



d^pi ( Hp 
Hp 



A^ = ie 2 e 2 N c J -^{_J^_Tr[ 75 (- i> 2 + m) Ap'ia + m) 7 ^i + m)] 

Tr[ 75 (- i> 2 + m^ifab + m) ^1 + m)]}, (2.17) 



where 



Pia=Pi~q, Pib = q~P2, 

N\ = p 2 — m 2 + ie, N2 = p\ — m 2 + ie, 

N ia = PL ~m 2 + ie, N ib = p 2 b -m 2 + ie, (2.18) 
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and e q is the electric charge of quark: e q = 2/3 for c quark and e q = —1/3 for b quark. The first 
and second terms in Eq. (|2.17p come from diagrams Fig. [2] (a) and (b), respectively. 

For the calculation of the form factor Fp^(q 2 ), it is convenient to choose the purely transverse 
frame q + = 0, i.e., q 2 = —q\ < 0. The advantage of this choice is that there is no the so-called Z- 
diagram contributions. The price is that only the form factor at space-like regions can be calculated 
directly. The values at the time-like momentum transfer q 2 > regions are obtained by analytic 
continuation. In this study, the continuation is not necessary because we only need the form factors 
at q 2 = for the P — > 77 and V — > P7 processes. 

At first, we discuss the calculation of Fig. [2Ja). The factors N\, N2 and iVj a produce three 
singularities in the p± complex plane: one lies in the upper plane; the other two in the lower 
plane. By closing the contour in the upper p± complex plane, the momentum integral can be 
easily calculated since there is only one singularity in the plane. This corresponds to putting the 
antiquark on the mass-shell. Given this restriction, the momentum P2 — > P2 with y\ — m 2 = 0, 
and pi = P — p2 ■ The on-shell restriction and the requirement of covariance lead to the following 
replacements: 

Nl -> N 1 =x 1 {M 2 -M^), 

N ia -► N ia = x 2 q 2 -x x M 2 + 2p ± -q ± , 

N 2 -> N 2 = N 1 + (1- 2xi)M 2 = x 2 M 2 - xiM 2 , 
d 4 pi Hp . f dx 2 d 2 p ± hp 

- m / /o-U AT" ( 2 - 19 ) 



(2vr) 4 N^Nia J (2vr) 4 x 2 N 1 N i( 

For Fig. |2](b) , the contour is closed in the lower p\ complex plane. It corresponds to putting 
the quark on the mass-shell and the momentum p\ — ► p\ with p 2 — m 2 = 0. In this case, we need 
to do the following replacements 

iV 2 -» N 2 = x 2 (M 2 -M 2 ), 

Nib -> N ib = x\q 2 - x 2 Ml - 2p L ■ q±, 

JVi -> iVi = xiM 2 - x 2 M 2 , 

d A pi Hp . f dx 2 d 2 p ± hp 

- ,7r / ^ m TT^T- ( 2 - 20 ) 



(2^) 4 NxN 2 Nib J (2^) 4 xiN 2 Nib 

From Eqs. (12.190 and (I2.20p . we see that Nn is obtained from Ni a by the exchange of x\ <-> X2 and 
the change of the sign of p± . 

After the above treatments, the transition amplitude of P — > 77* is obtained as 



pvo* f dx 2 d 2 P± N c e\m h P 
rP q x e J 



" ' 4tt3 Xl x 2 (M 2 -M 2 ) 

1 1 

+ 



-x 2 q 2 + XiM$ - 2p± ■ q± ' -x\q 2 + x 2 M^ + 2p± ■ q± 
Thus, the final formulae for the form factor Fpr 1 (q 2 ) is 

e 2 V2N~ c /■,„,, \ m 



(2.21) 



2 e Q ^J2N c f „ 
F Pi\<l ) = g7r 3 / dx 2 d p±(f)p{X2,P± 



m 2 + p\ 



+ 



x\Ml - x 2 q 2 - 2p ± ■ q ± x 2 M$ - x x q 2 + 2p ± ■ q± 



(2.22) 
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and Fp 7 (0) is 

e 2 \/2N r ml 

«•»«>) = - v^ / wVfrfe!i, ^sw (2 - 23> 

By comparing Eq. (|2.23|) with Eq. (|2.12|) , they share some similarities except the propagators N ia ^ 
(and a trivial factor e?). This point will become clearer when we discuss the relation between the 
light-front QCD and pQCD method. 

The decay rate for P — ► 77 is obtained from the transition form factors by 

r(P - 77) = ^(47ra) 2 |Fp7(0)| 2 - (2.24) 



D. V -» P7 



Similar to the analysis of P — > 77, we also consider a more general process of V — > P7* where 
the final photon is off-shell. The V — > P7* transition is parameterized in term of a vector current 
form factor V(q 2 ) by 

r„ = iee fluaf3 e"q a P (3 V(q 2 ). (2.25) 

where is the amplitude of V —> P7* process. P (e) is the momentum (polarization vector) 
of the initial vector meson, P denotes the momentum of the final pseudoscalar meson, and the 
momentum transfer q = P — P . To the lowest order approximation, the V — * P7* transition is 



7*(?) 




FIG. 3: Feynman diagram for V — > P7* process where P in the parenthesis denotes the momentum 
of initial meson and P denotes the momentum of final meson. 

depicted in Fig. [3j The amplitude is given by a one-loop momentum integral 

r — —iee N f f HvH ' p T I HyHp S b \e u (2 26) 
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where 



S% = Tr 



(jv - yy^iPl ~ P2)v) (- i>2 + m-hsW + m)7 At (^i + m) 



^ Wy 



— (Pi - P2^) (- p2 + rn)jfj,(- i>2 + m)7 5 (^i + m) , (2.27) 



and 



JV{ = p « - m 2 + i e; jV* = p ' 2 - m 2 + i e . (2.28) 



The first and second terms in Eq. (|2.26|) are arising from diagram (a) and (b) of Fig. [3J respectively. 
We have used the momentum relations: P = p\ + p 2 , P' = Pi +f>2> Q = P ~ P'i P2 = v'2 f° r diagram 
Fig. EKa); p\ = p' x for diagram Fig. (3](b). It is easy to find that S 1 " = S v . 

The momentum integral of Eq. (|2.26j) are performed analogous to the case of P — > 77*. The 
contour integrals are closed in the upper p~[ half-plane for the first term in Eq. (12.260 which 
corresponds to putting antiquark on the mass-shell; and in the lower half-plane for the second 
term which corresponds to putting quark on the mass-shell. For the first term, it leads to the 
replacements 

;vf - jvW = x 1 (mW 3 -mW 2 ) > 

d 4 pi H V H' P _^ _. f dx 2 d 2 p± hyh'p 

(2tt) 4 N!N 2 N{ ~* m J (2vr) 4 x 2 N 1 N[ 1 ' ' 

In order to preserve the covariance of the decay amplitude, we also need the replacements 

-,2 



P? - *iP a ~ 5°^, Vtvl - -j* (pi + ^0!) . (2.30) 



The similar treatments can be done for the second term. After using the above replacements and 
Eq. (|2.6p . we obtain the formulae for the form factor V(q 2 ) as 

2x _ e q [j A 2 4>v(x2,P±)4>'p(x2,P±) j 2 ( 2 (p± ■ q ± ) 2 ^ ] 



The rate for V — > Pj is 



III. NON-RELATIVISTIC APPROXIMATION AND PERTURB ATIVE QCD 

It is well-known that the system of the heavy quarkonium can be treated non-relativistically [l| ■ 
A relativistic invariant theory, light-front QCD in our case, should reproduce the previous results 
in the non-relativistic approximations. Here, we will explore the non-relativistic approximations of 
the light-front QCD. It is similar to the studies of heavy quark limit for heavy meson 0] within the 
light-front approach. In addition to it, the light-front QCD is related to perturbative QCD at the 
large momentum transfers, such as in P —* 77 process. Both of them show the different aspects of 
light-front QCD. 
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At first, we discuss the non-relativistic approximations of the light-front approach. In the rest 
frame of the heavy quarkonium, the momenta of quark and antiquark are dominated by their rest 
mass m 3> Aqcd (Aqcd is the hadronic scale). The momentum fractions x±,x 2 are peaked around 
|, and X2 — \ is of order Aqcd/"1- In NRQCD, the velocity of heavy quark is chosen as the 
expansion parameter. Neglecting the terms suppressed by 1/m, the invariant mass Mq and p z can 
be approximated as 

M ^ 2m M, Vz = (a?2 - \)M ~ Aqcd- (3.1) 

Compared with m, we have neglected the transverse momentum p± because it is of the order 
of Aqcd- Thus, the magnitude of the relative momentum p will be much smaller than m, i.e., 
\p\ ~ Aqcdj which constitutes the basis of the non-relativistic treatment. 

Under the non-relativistic approximations, the dependence of the hadron wave function on x 2 
is replaced by its dependence on p z since Mq = M is a constant. In this way, the hadron wave 
function will depend on the relative momentum p only, in other words, it can be represented by 
ip(p). The relation between the non-relativistic function ip(p) and the relativistic one (fr(x2,p±) can 
be established as follows. From Eq. (12.2|) . we obtain 

dp z M dx 2 , d 3 p = M dx 2 d 2 p ± , (3.2) 

As uausl, the function of ip(p) is normalized as 

d 3 p 



Comparing Eqs. p.2p and (|3.3p with Eq. (|2.8p . it is straightforward to derive a relation 

<P(x 2 ,p±) = V2Mi>(p). (3.4) 

Note that the above relation is valid within the non-relativistic approximation and is not correct 
in the general case. 

The hadron wave function in the coordinate space ^(f) is obtained by using the Fourier trans- 
formation 

d 3 p 

At the origin r = 0, ^(O) = / >^? a ip(p) is an important parameter which gives the magnitude 
of quark-antiquark coupling to the quarkonium. In the non-relativistic approximations, one can 
safely neglect p compared to m. For example, 



m = 1 *w ( 3 - 5 ) 



1 



m 2 + p\ — > m, X2 — — — > 0. (3-6) 



After these approximations, we can rewrite the decay constants Eqs. (|2.12p and (|2.14|) as 



/P .V^#, fv^VK^m. (3.7) 



Thus, 



f 2 p _ My |tfp(0)|' 



fy M P |^(0)| 2 ' 



(3i 
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This is just the so-called Van Royen-Weisskopf formula 26f] . Since the differences between vector 
and pseudoscalar vectors arise from the higher order in 1/m, these differences vanish in the limit 
m — ► oo, thus My = Mp = 2m, *&y(p) = *&p(p). The ratio of decay constants is equal to 1 in the 
limit. 

For the form factor V(0) of V — > -P7 process, Eq. (|2.3ip can be reduced to 

. f d 3 p 2^Mp-My^ v {p)^ P {p) 

Similarly, in the non-relativistic limit m — > 00, the form factor V(0) can be further written in a 
simple form as 

V(0) = 2e q /m. (3.10) 

Thus V(0) is a constant, independent of ^(0) because of the normalization condition of ^>(p). The 
physical picture is: the heavy quark and antiquark in the initial and final quarkonium are in the 
same momentum configuration at q 2 = point. It is analogous to the meson system with a single 
heavy quark that the Isgur-Wise function is normalized to 1 at the zero-recoil point in the infinite 
heavy quark mass limit. From Eqs. (|2.32p and (|3.10j) . the rate for the V — * Pj process is reduced 
into 

1 fi k 3 

T( y ^ p 7) = -ae\^. (3.11) 
where fc~ = (My - Mp)/2M V is the energy of the photon. This is the leading order result of Eq. 



121. 



(37) in 

Next, we discuss that the pQCD is applicable in P — > 77 process. In the rest frame of the 
heavy quarkonium, the total energy is 2m 3> Aqcd- Each final photon contains high energy of 
m and moves in the opposite light-front direction. When the high energy photon hits on one 
nearly rest constituent of the quarkonium, it causes a large virtuality of the order of m 2 . In 
particular, the virtuality of the internal quark is about 2m 2 from Eqs. (|2.19j) and (|2.20p . The 
transverse momentum in the propagator of the virtual quark can be neglected. Up to leading order 
in Aqcd/'' 71 ) the transition form factor i ? p 7 (0) is represented by 

dx2d 2 pj 
~(27r) s 

oc / dx 2 $(x 2 )T H (x 2 ). (3.12) 



F P7 (0) = e 2 q ^2N c J ( ^^- L 0(x 2 ,P±) T H (x 2 ) 



where &(x 2 ) is the hadron distribution amplitude obtained from wave function by integral over the 
transverse momentum, and Th(x 2 ) is the hard scattering kernel from the subprocess of qq — ► 77. 

The hard scattering kernel depends on momentum fraction x 2 when the loop corrections are 
taken into account. But, at tree level, the hard scattering kernel is 

T H = -L (3.13) 

It is not only independent of transverse momentum p±_ but also of longitudinal fraction x 2 . We 
thus have a further result 

F P7 (0) = e 2 ^. (3.14) 
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This equation means that the form factor i*p 7 (0) is proportional to the decay constant fp in 
leading order Aqcd/^ and leading order of strong coupling constant a s . After combining Eqs. 
(|2.15p . (|2.24p . (|3,7p . and (|3. 14[) . we finally obtain the decay rates for processes of V — > e + e~ and 
P — > 77 as 

r(y -> e + e -) = 1 4n c Wc v ^— 2 



3 My 



r(P^77) = 16We*^|^. (3.1! 
These results are the same as ones in Table III in the non-relativistic quark-potential model [hJ 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In order to obtain the numerical results, the crucial thing is to determine the momentum dis- 
tribution amplitude (p(x2,p±)- One wave function that has been often used in the literature for 
mesons is the Gaussian-type 

with N = 4(vr/^ 2 ) 3 / 4 and 

^ = * . (4.2) 

d,X2 X\XiMq 

The required input parameters include quark mass: m c for c quark and m& for b quark; a hadronic 
scale parameter (5 for rj c ^) and J/ip(Y). The quark mass entered into our analysis is the constituent 
mass. For light quarks (u and d), the constituent mass which several hundred MeV, is quite bigger 
than the current one which is only several MeV obtained from the chiral perturbation theory. While 



for the heavy quarks, the difference between them is small. From PDG [27j], the current masses are 
1 GeV < m c < 1.4 GeV and 4 GeV < mb < 4.5 GeV in the MS renormalization scheme. For our 
purpose, we will choose heavy quark constituent masses as 

m c = 1.2 GeV, m b = 4.3 GeV. (4.3) 



Our choices are smaller than the parameters given in [8j , but they are consistent within the error of 
one Aqcd- For the meson mass, M Vc = 2.980 GeV, Mj^, = 3.097 GeV and My = 9.460 GeV Q. 
The mass of % is still unknown and it is parameterized as Am = M-f — M Vb . From the references 
in [ijj], the range of Am is Am = 30 — 150 MeV. 

After fixing the quark and meson masses, the remained thing is to determine the parameters (5. 
For the vector meson, fly is extracted from the decay constant fy which is obtained directly from 
the process V — ► e + e~ by Eq. (|2.15|) . For the pseudoscalar meson r/ c , (3 Vc is extracted from the 
decay constant f Vc which is obtained from the process B — > r] c K. 

For the cc charmonium system, there are some experiment data which provides a place to test 
the applicability of the Gaussian-type wave function to the heavy quarkonium. From J/ij) — ► e + e~, 
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we obtain fj/^ = 416 ±6 MeV, and extract = 0.639 ±0.006 GeV. Prom B — > rj c K, one obtains 



= 335 ± 75 MeV [28], and we extract (3 Vc = 0.652 



+0.165 
-0.143 



GeV. It is apparent that the dominant 



errors in the following calculations will be derived from the uncertainty of f^ c . By using the above 
parameters, we give the numerical results for r\ c — ► 77: Br{jj c — > 77) = (1.78 ~ 3.05) x 10 -4 and for 
J/t/j — > ^7: Br(J/ip — ► ^7) = (2.38 ~ 2.84) x 10 -2 . The experimental data are Br(rj c — ► 77) = 
(2.8 ± 0.9) x 10 -4 and Br(J/^ — > 7^7) = (1.3 ± 0.4) x 10~ 2 . Obviously the former fits experiment 
very well but the latter does not. This inconsistency still exists even we adjust the quark mass m c 
in the range 1 ~ 1.4 GeV. 



One may consider a power law wave function similar to the one employed in Ref. 29J to fit the 
data, however, the Gaussian-type wave function has been used widely in the phenomenal analyses 
which related to meson. Thus we modify the Gaussian-type wave function by just multiplying a 
factor (x\x 2 ) n 

„2\ 



<t>(x 2 ,p±) = N( Xl x 2 ) n \ 



dp z 



exp 



P± +P Z 



(4.4) 



, dx 2 \ 2(3 2 , 

The curve which x 2 is peaked around | will be sharped or dulled if n > or n < 0, respectively. 
In the non-relativistic limit, Eq. (|3.6[) reveals that the curve is near to a delta function 5(x 2 — |). 
Therefore the case of n > seems suitable for the heavy quarkonium. In fact, if n = 5 and m c = 1.2 
GeV, we can extract j3j /l/} = 0.786 ± 0.008 GeV and f3 Vc = 0.807lg|[f GeV. The numerical results 
Br(r] c -> 77) = (1.56 ~ 2.06) x 10" 4 and Br(Jj^ -> ^7) = (1.62 ~ 2.41) x 10 -2 are both consistent 
with the experimental data. Thus there is a deduction that, for heavy quarkonium, the momentum 
fraction x 2 is more centered on i than one is in the Gaussian-type wave function. We show the 
x-dependent behaviors of these two types of wave functions in Fig. [4] and the numerical results in 
Table E 




FIG. 4: The x-dependent behaviors of <fr (dash line) and <j> (solid line, n=5) at p\ = 0.1 GeV 2 



For the bb bottomonium system, the experimental data are relatively less. From T(1<S) — > e + e , 
we obtain / T = 708 ± 8 MeV, then extract fa = 1.323 ± 0.010 GeV and (3 r = 1.463 ± 0.012 GeV. 
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TABLE I: The comparisons between the experimental data and theory predictions for charmonium 
decays. 





Br(r) c -> 77) 


Br{J/ip -> r/ c 7) 


experiment data 


(2.8 ±0.9) x 1(T 4 


(1.3 ±0.4)% 


this work (<fi) 


(1.78 ~ 3.05) x 10~ 4 


(2.38 ~ 2.84)% 


this work (0) 


(1.56 ~ 2.06) x 10~ 4 


(1.62 ~ 2.41)% 



However, rjb meson hasn't been observed in experiment. It is impossible to determine the decay 
constant /„, from the experiment. As had been discussed, the relation f„ b = ff is hold in the 
non-relativistic limit. Since the corrections to this relation are suppressed by Aqcd/^&j it may be 
reasonable to use it to determine the parameter /3„ (/?„ ). Therefore we obtain (3^ = 1.433 ± 0.014 
and Prj b = 1.607 ± 0.018 GeV. For obtaining the decay widths T(r)b — ► 77) and T(T — ► 7/yy), we 
must be aware of the value of Am. However, the sensitivities of these two decay widths to Am 
are quite different. On the one hand, r(r/& — ► 77) is insensitive to Am because M % 3> Am (see 
Eq. (|2.24p ). On the other hand, T(T — ► 7^7) is very sensitive to Am because it is proportional to 
(Am) 3 (see Eq. (12.32P ). Thus here we list the values of T(rjb — > 77) for Am = 0.09 ± 0.06 GeV and 
r(T — > Tjbj) for Am = 0.09 GeV in Table [Hi The dependences of T — > 77^7 on Am are also shown 
in Fig. [5j 



TABLE II: The comparisons among the several theory predictions for bottomonium decays. 

rjjjb ^ 77) (eV) r(T -> wy) (eV) 

this work (0) 453 ± 17 33.2 ±0.1 

this work ($) 422 ± 15 31.5 ±0.1 

used in [14] 557 ± 85 

NRQCD [30]O(a s ) 460 
potential model [31j 466 ± 101 



For the numerical results, some comments are in orders: 
(1) The decay constant for rj c is f Vc = 335 ± 75 MeV [28], and we obtain 

2 

i w 0.65 ±0.31. (4.5) 



fv 



\ fj/i> J 

The difference between the pseudoscalar and vector meson in light-front approach comes from 
power suppressed terms: the transverse momentum p±, X2~ \ and the wave function 4>(x,p±). The 
deviation of the results from 1 shows that Aqcd/^c ~ 30% corrections cannot be neglected. 

(2) For the Ml transition J/ip — > r/ c 7, the leading order prediction from Eq. (I3.1ip for the 
branching ratio is 4.2% which is about a factor of 3 larger than the experimental data. This means 
that the next-to-leading order Aqcd / Tn c corrections are so substantial that they must be included 
in the calculations. 
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FIG. 5: The dependences of T(T — ► 7J&7) on Am = M~f — M, 



(3) For the decay width r(r?b — ► 77), it is insensitive to the variations of Am but proportional 
to f%. Thus an adjustment of f Vb by 10% will correspond to a variation of the theory prediction 
by about 20%. So far, (f Vb /fr) 2 = 1 is our assumption and 0.99 ± 0.04 in [3j and 1.16 ± 0.06 in 



33J. After considering these uncertainties, our prediction may be consistent with previous results 
listed in 0, HE I- 



V. CONCLUSIONS 



In this article we have studied the decay constants, two-photon annihilation P — > 77 and 
magnetic dipole transition V — * P7 processes for the ground-state heavy quarkonium within the 
covariant light-front approach. The phenomenological parameters and wave functions are deter- 
mined from the experiment. The predictions agree with the measured data within the theoretical 
and experimental errors. The quark mass we use is very close to the current mass which is different 
from the choices in non-relativistic quark model. The difference between the J/ip and r\ c decay 
constants and the study in J/ip — ► r] c j both show that the power corrections from wave functions 
and transverse momentum effects are important. In order to make a better fit to the experimental 
data, we adjust the longitudinal momentum fraction of the wave function to center around 1/2 
further. We also give a numerical prediction for r]b — * 77 and T — * rjb'y. The branching ratio for 
~f ~~ > %7 is too small to be observed. 775 — ► 77 may be a good process to determine rib and its 
mass. The QCD corrections are neglected in this study, including them will slightly change the 
wave function inputs but does not change our conclusions. 

The light-front approach shows different aspects of QCD. Under the non-relativistic approxima- 
tions, the light-front approach reproduces the results in the non-relativistic quark-potential model. 
For P — > 77 where two final photons are on the opposite light-front, the process is perturbative 
dominated and the light-front approach reduces to the model-independent pQCD. The light-front 
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method unifies the perturbative and non-perturbative QCD into the same framework. 

We have considered s-wave heavy quarkonium in the light-front approach only, the applica- 
tions to other quantities and higher resonances are in progress. One interesting thing may be 
to explore the light-front approach in NRQCD (or pNRQCD). This will provide an alternative 
non-perturbative method to calculate the hadronic matrix elements defined in NRQCD. 
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